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Abstract 

In the presence of a sharp corner in the boundary of the entanglement region, the en¬ 
tanglement entropy (EE) and Renyi entropies for 3d CETs have a logarithmic term whose 
coefficient, the corner function, is scheme-independent. In the limit where the corner becomes 
smooth, the corner function vanishes quadratically with coefficient cr for the EE and for 
the Renyi entropies. For a free real scalar and a free Dirac fermion, we evaluate analytically 
the integral expressions of Casini, Huerta, and Leitao to derive exact results for a and (t„ for 
all n = 2, 3,.... The results for a agree with a recent universality conjecture of Bueno, Myers, 
and Witczak-Krempa that aIC t = 7r^/24 in all 3d CETs, where Ct is the central charge. For 
the Renyi entropies, the ratios (JuICt do not indicate similar universality. However, in the 
limit n ^ oo, the asymptotic values satisfy a simple relationship and equal l/(47r^) times the 
asymptotic values of the free energy of free scalars/fermions on the n-covered 3-sphere. 



1 Introduction and Results 


For a 3d conformal field theory (CFT) in the ground state, the entanglement entropy S for a 
region whose boundary has a sharp corner with angle 9 can be written as 

5 = -a(0)ln(^)+0(1). (1.1) 

Here L is a length scale associated with the size of the entangling region, e is a short distance 
cutoff, and H is a non-universal constant. The corner contribution to the entanglement entropy 
is the scheme-independent positive function a{9) of the opening angle 9 [1, 2, 3]. Since the 
entanglement entropy of the region equals that of the complement region, the corner contribution 
satisfies a(27r — 9) = a{9). If the curve bounding the entangling region is smooth, the logarithmic 
term is absent, hence a{9) must vanish in the limit 0 —>■ vr and it does so quadratically as 

a{9) = a {9 — tt)'^ + ... for 0 —)■ vr. (1-2) 


The value of the corner coefficient a depends on the theory. 

For the theory of a free real scalar or a Dirac fermion, Casini, Huerta, and Leitao [4, 5, 2] 
derived expressions that give a{9) implicitly in terms of some rather involved integrals. In the 
limit, 9 ^ TT one can extract double-integral expressions for the corner coefficient a in (1.2). These 
integrals have been evaluated numerically [5, 6] and the results indicate that the exact values are 
[ 6 ] 

(1.3) 

for the free boson and free fermion, respectively. 

Bueno, Myers, and Witczak-Krempa [6] conjectured that the ratio of the coefficient a in (1.2) 
to the central charge Ct is universal in 3d CFTs and that it takes the value 


= - and = - 

256 128 


conjecture [6]: 


a 


IT 


Ct 24 ■ 


(1.4) 


The conjecture (1.4) has passed non-trivial holographic tests for gravity models with a family 
of higher derivative corrections [6, 7]. The central charge Ct is defined as the coefficient of the 
vacuum 2-point function of the stress tensor (see eq. (3) in [6]). For free bosons and fermions, 
Osborn and Petkou [8] found that C^^ = 3/(327r^) and C^'^ = 3/(IGtt^) in 3d. So with the values 
(1.3), the ratio a (Ct is indeed 7r^/24 for both free bosons and fermions. 

In this paper, we evaluate analytically the integral expressions [6] of Casini, Huerta, and Leitao 
[4, 5, 2] for and and prove that their exact values are indeed those in (1.3). This verifies 
the universality conjecture (1.4) for the case of free bosons and fermions. One way of viewing 
the conjecture is simply as the statement that the corner coefficient cr in (1.2) does not contain 
independent information about the CFT, but is fixed in terms of the central charge Ct- 

Turning to the Renyi entropies Sn^ one can define a similar corner contribution 0^(0) which 
in the smooth limit 0 —)■ vr goes to zero as a„(0) = cr„ (0 — vr)^ + ... for n = 2, 3,4,.... (The 
n —)• 1 limit of the Renyi entropy is the entanglement entropy.) It is not known if g^ICt has any 
universal properties. 
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We calculate a-n analytically for the free boson and free fermion using integral expressions for 
an derived in [4, 5, 2].^ For the free scalar we find 



(B) ^ A:(re - A:)(re - 2A:) tan (^) 

4^ 247rre^(re —1) 

k=i ^ ’ 

(1.5) 

Note that when re is even, the contribution from k = re/2 must be taken carefully using 
limfc^„/ 2 (re - ‘Ik) tan (^) = 2re/7r. 

The result for the free fermion is 


(F) fc(re2 - 4A;2) tan (^) 

(1.6) 

^ 247rre^(re—1) ’ 

fc=-(n-l)/2 ^ ' 


where sum is to be taken in integer steps from — to 

For low values of re, the finite sums of the trigonometric functions in (1.5) and (1.6) simplify 
quite nicely. The results for first nine values of an are 


re 


a 


[B) 


2 

3 

4 

5 

6 

7 

8 

9 

10 


48-712 


0.00211086 


108\/377 

0.00170163 

8+377 . 

1152772 ~ 

0.00153255 

\j2'Q-2s/l 

100077 

Ri 0.00144219 

81+34^/377 

19440772 

R7 0.00138643 

2cot(^)+5 cot( j|)+5tan(y) 

411677 

32+9^(l+2V2) ^ 


0.00134874 


107527r2 

27^3+10 cot( )+28 tan( f )+35 tan( ^) 
349927 


0.00130116 


a 


(F) 


1 

647r 


0.00497359 


5 

216\/37r 

1+6^2 

7687r 


77S 0.00425408 
Ri 0.00393133 


-v/425+58V5 

20007r 


0.00374840 


261+20v4! 

259207r 


0.00363061 


13cot(^)+22cot(|^)+15 tan(y) 
82327 


0.00354841 


= «.«0S4S777 

135^3+68 cot(^)+77tan(f )+130tan{^) 
6998477 


0.00344118 


5+300\/5+4\/425+58\/5 

7200077 


0.00340427 


In the case of the scalar, the exact re = 2,3 results were guessed by the authors of [6] based on 
their high precision numerical evaluation of the integrals. 

Since the ratios of the central charges of free fermions and bosons differ only by a factor of 
2, universality of the ratio anjCT would require that a^ /an obeys some simple, possibly re¬ 
dependent, relation. Based on our results above, there is no hint of such a simple relationship. Of 
course to fully exclude this, one would need values of an for other 3d CFTs. 

As a function of re, the Renyi corner coefficient cr„ decreases monotonically, as shown on the left 

We are grateful to Horacio Casini for sharing with us the integral expression for cr„ . 
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crn/o'„-(1 + 1/n) 



Figure 1: Left: Plot showing that cr„ decreases monotonically from the entanglement entropy value included 

for n = 1 to the asymptotic value ctoo for free scalars (blue) and free Dirac fermions (maize squares). The 

asymptotic values « 0.0011569 (black) and « 0.00308507 (gray) are indicated as 

horizontal lines. Right: The plot illustrates our numerical fit (t„ = croo(l + ^ + ^ + ^ + -- -)’ which 
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we find 6 i = 62 = 63 = 1 for the free scalar, and 61 = 1 and 62 = ^3 = 1 — 12 ^( 3 ) ~ 0.31578 for the free 
fermion; the solid curves are ^ ^ for those respective values of 62 and 63 . 


in figure 1. When n is large, an asymptotes to a constant value, which we calculate analytically: 


O', 


{B) 


3C(3) 

327r4 


and 


O', 


{F) 


m. 


(1.7) 


The appearance of the Riemann zeta-function is intriguing since C(3) ~ 1.20206 also shows up 
in the free energies and Renyi entropies for free scalars/fermions on a 3-sphere, as shown by 
Klebanov, Pufu, Sachdev, and Safdi [9]. Specifically, the free energy of a free real scalar or free 
fermion on an n-covered 3-sphere behaves as^ J-n —t nJ-^o for n —?• 00 with 


and = (1.8) 

Thus, for both free scalars and fermions we have 


^ 00 


J_ jr{B/F) 

47r2 “ • 


(1.9) 


For finite n, there is no apparent relation between fFn and an, however there are some similarities 
in the subleading large-n behaviors, as we discuss in section 4. The plot on the right in figure 1 
shows the large-n behaviors of the Renyi corner coefficients an- A priori it is not clear if there is 
any relation at large n between an and Fn, but it would be curious to test (1.9) in other examples. 

The remainder of the paper details the derivations of the results summarized above. In section 
2, we derive the results (1.3) for the entanglement entropy corner coefficient a. We then evaluate 
the Renyi entropy corner coefficients an in section 3. In section 4, we discuss the asymptotic 
behavior at large n. 

^The authors of [9] work with a complex scalar, so the free energy there is twice that of a real scalar. 
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2 Evaluation of the EE integrals 


In this section we describe the procedure for analytically evaluating the integrals for the coefficients 
and of the entanglement entropy. Our starting point is the integrals [4, 5, 2] presented 
in equations (B1)-(B3) of [6]. After a change of integration variable from m to /r = \/4m^ — 1, 
the integrals take the form 





db (1 — a) 


IT 


cosh^ (vrb) 


db 


fi^Ha (1 — a) — 


flF TT 
dvr sinh^ (tt^) ’ 


( 2 . 1 ) 


where a = 1/2 + ib for the scalar and a = ib for the fermion. The functions H and F are defined 
as 




1 


Ibvra (a — 1) ’ 


F = -^ 
F2 


with^ 


Fi = AnchHa {1 — a) (2a — 1)^ +—-c/i^ (/x^ + l) 


ch 


F2 = 


(2o - 1)^ + 


1 2 


2 (2a — 1) 

The functions h, c, Xi, X 2 , and T are defined as follows: 
2 (lU^ + (2a — 1)^) sin^ (ira) 


h = 


c = 


+ 1) (cos (27ra) + cosh (vr/i)) ’ 

2^“7ra (1 — a) sec ^vra + T — a + 


Ai = - 


v/)?TT (r (2 - a))" r(-i + a + f) 
r(— a) TTsinh (^) + icosh (^) 


i+a+f i + a-f 


( 2 . 2 ) 


(2.3) 


A 2 

T 


22a+i^r(a + 1) — a + T^^ — a — (cos (27ra) + cosh (vr/x)) 

(Ai with a replaced by (1 — a)) , 


^\/h[{h + 1) (/x2 + 1 ) - 4a (1 - a)] 


(2.4) 


Here denotes the digamma function, 't/{z) = ^ Inr(z). 

Our first line of attack involves calculating the quantities cXi/h and X 2 /C that appear in H 
in (2.2). Beyond the immediate cancellations that occur in these ratios, one can perform further 
simplifications using identities involving gamma functions. Namely, one can use the recurrence 
relation 

r(i + z) = zr(2:) (2.5) 

^We simplified the expression for F\ in [6] by writing it in terms of H. 
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and the reflection relation 


r(i-z)r(z) 


vr 

sin (vrz) 


Surprisingly, all the gamma functions cancel after a series of such substitutions, giving 


c 

h 


Xi 


1 

c 


^2 


1 CSC (vra) 
IdTTfi {a — 1) a 


vr sinh (^) + i cosh 


\/?T 1 CSC (tto) 
IGttiU (a — 1) a 


vrsinh (^) +icosh (^) 



+a- f)) 

f)) 


( 2 . 6 ) 


(2.7) 


It is suggestive that the pre-factors and the form of these two results are the same. We then 
proceed by adding them together as in (2.2). The linear combination of digamma functions that 
appears in the result can be simplified using properties easily derived from (2.5) and (2.6). In the 
form that is useful for our purpose, these identities are 

f) = f) + 1 _ ^ ^ 

2 ^ ^ 2 

and 

+ a ± ~ ® ) = vrtan (vra ± • 

Then the combination of Xi and X 2 that appears in H simplifies to 


CSC (vra) cosh (^) \ 

(l-2a)2 + ^2 j 

The last ingredient we need to construct H in (2.2) is T. Using (2.4), it is 


I^ _ y^u2~+I / 7rsin(7ra) sinh (^) 

h c 47ra(I — a) \ fi [cos(27ra) + cosh(7r/r)] 


T = 




((1 — 2a)2 -|- ^ 2^2 cosh^ (^) 

{p? + 1) (cos(27ra) -t- cosh(7r^)) 


( 2 . 8 ) 


(2.9) 


Further simplifications of H depend on the nature of variable a, as we will see when we specialize 
to the cases of the free scalar and the free fermion. 


Free scalar. To proceed with the evaluation of the integral we set o = l/2 + i6 as prescribed 
for the free scalar. It is furthermore convenient to change integration variable b —)• 6/2. Using 
that both /i and 6 are positive, the integrand of simplifies dramatically and becomes 




f°° fj, [tt (/r2 — 52 ) sinh(7r/r) -|- 2/xcosh(7r6) — 2^cosh(7r/x)] 

/o ^ Jo 64[cosh(7r6) — cosh(7r/i)]2 


( 2 . 10 ) 


Next, we integrate by parts. Writing 


= - 
64 


poo PC 

/ d^i 
Jo Jo 


db 


1 /■“ , „ d ( (/^^ - 6^) 


+ 


62-/r2 


dfj, \ cosh (irb) — cosh (vr/r) j cosh (7r6) — cosh (tt^u) 


( 2 . 11 ) 
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we see that the boundary term vanishes and we get 


<7 


(B) _ 


1 


r+oo 


r*H-oo 




db 




( 2 . 12 ) 


256 J-ao cosh (irb) — cosh (vr/r) 

We have extended the limits of integration to facilitate the change of integration variables 

fi = X — y and b = x + y. (2-13) 

This separates the two integrations and reduces the expression to 




r-+0O 


dx 


sinh (ttx) 


1 


(2.14) 


This completes the derivation of the result (1.3) for the free scalar. 


Free fermion. With Fi given in terms of H as in (2.3), we have already done most of the leg-work 
needed to compute For the free fermion, we have to take a = ib and it is again convenient 

to change integration variable b —)• 6/2. After putting everything together, we have 


32, 


dy, 


y[7r {^y? — 6^ — l) sinh(7r^) — 2//cosh(7r6) — 2/rcosh(7r//)] 


'0 


[cosh(7r6) + cosh(7r/i)]^ 


(2.15) 


We can express the integrand as a total derivative plus remaining terms as 


= - 


1 1*00 pc 


db 


^ ( d -6^-1) 

dy I cosh(7r6) + cosh(7r/x) j ' cosh(7r6) + cosh(7r^) 


+ 


1 - /i2 + 62 


(2.16) 


As before, the boundary term vanishes and we are left with the expression (after extending the 
limits of integration) 


u(^) = 


1 


/ -l-oo 

-oo 


dy db 


1 — y? Fb"^ 


1 


r>-|-0O 


1 -k Axy 


cosh(7r6) + cosh(7r|u) 128 


dxdy- 

oo cosh(7rx) cosh(7ryj 


(2.17) 


In the last step, we changed integration variables using (2.13). Since x/ cosh(7rx) is odd, that part 
of the integral vanishes and the result is therefore simply 




1 

1^ 


^+00 \ 2 
/_oo cosh(7rx) 


1 


(2.18) 


Thus we have derived the result (1.3) for the free fermion. 


3 Renyi entropies 

We now proceed to calculate the corner coefficients an for the Renyi entropies. 

Free scalar. For the scalar field, the Renyi corner coefficient is given by the integral (B7) in [6]. 
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We change of the integration variable m to ^ Am? — 1 to write it as 


n-l 


= - 

^ 'T) 


k {n — k) f°° 2 rr 

^ 2n2 (n - 1) io ^ ’ 


(3.1) 


where is H in (2.2) with a replaced by k/n. With the simplified expression for H from 

section 2, we get 


r(S) = 


sm 


'TTk\ 
K n ) 


327rn2(n — 1) Jq 


dn 


jjL [(n — 2k)‘^ + iJ?n?~\ 7rsinh(7r/r) — 2iJ?v? [cos (^) + cosh(7r/i)] 


[cos (^) + cosh(7r//)]' 


(3.2) 


As before, we write the integrand as a total derivative plus the remaining terms: 


(iik\ 


dfj, 



^ fj.[{n — 2A:)2 -|- fj?r?'^ \ 

(n — 2k)‘^ + fj?n^ 


i^cos (^) + cosh(7r^) y 

cos (^) + cosh(7r/x) 


^ 327rn2(n- 1) Jq 

The boundary term vanishes and the expression simplifies to 

(n — 2k)‘^ + 


(3.3) 


= 


n—1 ■ 2 /Trfe' 

sm ' 


E 


(f) 




^ 327rn2(n - 1) do cos (^) + cosh(7r/r) ' 
The contribution of A: = n/2 is easy to calculate and is equal to 


(3.4) 


1 


d/i — 




1 


647r(n-l)7o ^ sinh^ (^) 487r2 (n - 1) 

For k 7 ^ n/2, there are contributions from two integrals: 

d/x _ 2tan“^ (tan (^)) _ 2 J 


(3.5) 


f 


= 

I" cos (^) + cosh(7r/x) 


TT sm 




sm 




^ A: < n/2 
^-1, A:>n/2 


(3.6) 


and 


r(2) _ 

^n\k 


f 


/x^ d/x 


[ / 2ik7r \ / _ 2ik7T \ -i 

Li 3 ( - e " j - Li 3 ( - e - jj 


cos (^) + cosh(7r/x) 

Uog (e^)[7r2 + log2 (e^)] 


TT"’ sm 


in(^) 


(3.7) 


k(rd-ikJ) 

3n3 


A: < n/2 


37r3sin(^) sin(^) [ _ (n-fc)(n-^fc)(3n^ ^ k>n/2. 

Above, we manipulated the tri-logarithm Li 3 using the polylog identity 

Li 3 (z) -Li 3 (w^) = -^log^(-z) - y log(- 2 :) , (3.8) 

which holds for z ^ ]0,1[. 

Combining the results (3.6) and (3.7), we find that the result is the same for 1 < A: < n/2 and 
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n /2 < k < n, namely 


(3.9) 


dn 


(n — 2k)‘^ + 


= („-2km+n^& = Sk(n-k)(n-2k) 


Jo cos (^) + cosh( 7 r^) 3n sin(^) 

/ D\ 

Thus, having evaluated the integral in (3.4), we can write an as the finite sum 

1 


n—l 


= 


—- 5 —-- kin — k)(n — 2k) tan(^) . 

247rn3 (n - 1 ^ ^ ^ ^ 

^ ’ k=l 


(3.10) 


Note that taking the limit k —>■ n/2 as described below (1.5), the summand evaluates precisely to 
the special case (3.5). The expression (3.10) is the result for the Renyi corner coefficient presented 
in (1.5), so this completes our evaluation for the free scalar. 


Free fermion. For the fermion field, the Renyi corner coefficient is given by the integral 

/■“ r //FI 

'y ldn a{l-a)^i^H- — 

'0 I ^'^la=k/n 


= - 
'n 


n — 


I E 

k>0 


(3.11) 


where the sum is over k from 1/2 (n even) or 1 (n odd) in integer steps to ^(n — 1). Substituting 
the expressions for H and F obtained earlier gives 




2 ^" sin^ (—) f°° [cos (^) + cosh(7r^)] - /i - 1^ 7rsinh(7r//) 


E 

fe >0 


87r(n - 1) 7o 


I 




[cos (^) +cosh(7r/r)]' 


(3.12) 


We then use integration by parts to simplify the integral 


CJ 


F) 


|(n-l) 

E 


k>0 


(w) 

87r(n — 1) 




f k 


) 


dpi 

i cos 

^)+cosh(7rM) J 

cos (^) + cosh(7r/x) 


(3.13) 


The boundary term integrates to zero and the expression simplifies to 


= 

n 


E 

k>0 


Sin 




dfi 


l-/i2_ 


4fc2 


87r(n — 1) Jo cos (^) + cosh(7r/x) 


(3.14) 


The result of the integral again involves a difference of two tri-logarithms and it can be simplified 
using equation (3.8). The result is even in /c —)• —k and we can write the final answer as 


CJ 


{F) 


247r n3 (n — 1) ^ ^ 

fc=—|(n— 1 ) 



(3.15) 


This is the answer we presented in (1.6). Values for low n were tabulated in section 1 for both 
and a^^. 


















4 Asymptotic behavior of the Renyi’s 


Let us now study the large n behavior of the Renyi entropy corner coefficients an- In particular, 
we evaluate analytically the value for the coefficients an in the limit where n —)• oo. This is done 
by introducing a new variable x = k/n and multiplying by nAx = 1. Then in the n —>■ oo limit, 
the sum becomes an integral and we have 


^ ^ dx x{x-l){2x- 1 ) tan (vrx) = , 

^ i /_i/2 • 


(4.1) 


These values turn out to be proportional to the asymptotic values of the Tn nToa calculated 
on the n-covered 3-sphere [9]; as noted in (1.9) we have . 

On the right in figure 1, we illustrated the asymptotic behavior of the corner coefficient which 
we find to be 

<T„ = <To„(l + ^ + h + l| + ...). (4.2) 

\ n n vr / 

Numerical fits show that bi, b 2 , and 63 are 1 for the free boson while 61 is 1 and 62 = ^3 ~ 0.31578 
in (4.2) for the free fermion. In fact, fitting up to 0(l/n^®), we find numerical evidence that 
= ^ 2 fc-i-i for both the scalar and fermion. This indicates that a factor of (n -|- l)/n can be 
factored out of the function in (4.2), so that 


O'n. — CTn 


re + 1 


re 


1 + 


62 ,^4 , be 


H—T H—^ + • • • 


re^ 


re’ 


re" 


(4.3) 


It is also interesting to study the ratios of the Renyi corner coefficients at large re: based on 
numerical fits in the range re = 100 to 2000 we find 






(F) 8 




re" 


(4.4) 


The value of the 1/re^-coefficient is guessed based on the numerics, 
function 


3 

8 


di d2 ds 
1 + — + ^ + ^ + ... 


re 


re^ 


re" 


Specifically, we fit to the 
(4.5) 


and find that di < 10“^®, |d 2 — 1 ^( 3 ) | < 10“^^, d^ < 10“^®, ^4 = —0.93871149..., ds < 10“^^ 

etc. The vanishing of the odd powers in (4.5) is consistent with (4.3). Note also that we can now 

identify the number 62 = ^3 ~ 0.31578 from the fit (4.2) of the free fermion Renyi entropy corner 

2 

coefficient at large re as 1 i2C(3); this is the value given in the caption of figure 1 . 

Taking the Hurwitz zeta-function expressions for 7^ ' from [9] and using (4.5) to perform 
a similar fit at large re in the range 30 to 300, we find 


7r(S) 

n 

'f n 



TT^ 1 

12C(3)re2 


+ 0.937106586^ O 



(4.6) 
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Again, the value of the 1/n^-coefficient is guessed based on the numerics which give di < 

\d 2 + ^ 2 ^( 3 ) I < 10“^'^, ds < 10“^^, d 4 = 0.937106586..., ds < 10“^^ etc. The behaviors of 
individually is, however, very different that that of the Renyi corner coefficients. We find that 
~ nT'i\\ + 0{^)) while J-f ^ ~ + 0{^)). 

It is not clear whether the similarities observed at large n between an and Fn have any 
significance or if it is a coincidence. Perhaps future investigations will clarify this. 
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